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ABSTRACT 

We present the first results of a Spitzer IRAC (Infrared Array Camera) wide-field survey of the Coma clus- 
ter The observations cover two fields; the first is a 0.733 deg^ region in the core of the cluster (Coma 1), the 
second a 0.555 deg^ off-center region located ~57' (1.7 Mpc) south-west from the core (Coma 3). The IRAC 
observations, although short 70-90 s exposures, are very sensitive; we detect ^29,200 sources at 3.6/im over 
the total ^1.3 deg^ survey area. We construct 3.6/im galaxy luminosity functions (LFs) for each field using 
selection functions based on spectroscopic redshifts. At the bright end, the LFs are well modeled by a tradi- 
tional Schechter function; (M5.6^„i,ai> = (-25.17, -L18) and (-24.69, -L30) for Coma 1 and Coma 3 
respectively. However, at the faint end (Ms.g^Tji > —20.5), there is a steep increase in the LF slope in both 
fields indicative of large numbers of red dwarf galaxies. The reality of this population is examined using SDSS 
optical counterparts with optical color filtering (g — r < 1.3). The steep increase can be modeled with a power- 
law function, with slopes of a2 = -2.18 (Coma 1) and 02 = -2.60 (Coma 3), the difference likely indicating a 
change in environmental effects between the two fields. A qualitative comparison with optical (B- and i?-band) 
LFs shows that we are likely to be observing a larger population of dwarf galaxies in the near-IR, which may 
be a low-surface-brightness (LSB) population that IRAC is particularly sensitive to, or a population too red to 
be detected in existing optical surveys down to i? ~ 20. 

Subject headings: galaxies: clusters: general — galaxies: clusters: individual (Coma) — galaxies: luminosity 
function, mass function 



1. INTRODUCTION 

The study of the luminosity function (LF) of galaxies in 
clusters is important to understand the formation and evolu- 
tion of galaxies. For example, a testable prediction of stan- 
dard cold dark matter (CDM) hierarchical clustering mod- 
els is a steep mass function resulting from the presence of a 
large number of faint dwarf galaxies, a fossil record of small 
dark matter halo formation in the early Universe (e.g. [White 



|& Rees|[T978l [Kauffmann eTaL][T993] l. The LF is deKried 
as the number density of galaxies per unit luminosity, and 
can be used to directly measure the abundances of various 
galaxy populations in different environments (e.g. Trentharn 
[& Tully,,2002) . In clusters, it can generally be described by 
a Schechter function (Schechter 1976| l, with the bright end 



(Mr ^ -18) dominated by giant elliptical galaxies, and the 
faint end (M^ > -18) dominated by red dwarf elliptical (dE) 
galaxies (e.g. Seeker & Harris'1996| Seeker et al. 1997'; Tren-J 



[fcam ^l998j Mobasher et al. 2003) . Observations of nearby 
clusters allow us to probe the faint end of the LF in detail, and 
there is now substantial evidence that clusters contain signifi- 



cantly higher numbers of low-luminosity dwarf g 


alaxies than 


are found in the field (e.g. 


Trentham et al.|2005 Popesso et 


|al.|2006HAdami et al..2006 





Although it seems that the most massive galaxies in clusters 
formed at early times (z > 2), it is not known whether dwarf 
galaxies also formed early-on within the cluster or whether 
they were accreted from the field at lower redshifts (z < 0.5). 
There is, however, growing evidence from infrared (IR) LF 
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Studies in support of the l atter scenario (e.g. de Propris et al. 
1999[ Muzzin et al.|2 005'). This is further supported by stud 



ies of the spatial distribution and radial velocities of dEs in 
the Virgo cluster ( [Conselice et al. 2001), where it was found 
that the cluster dwarf population probably originated as field 
spiral galaxies that were accreted an d transformed into d Es by 
the process of 'galaxy harassment' ( Moore et al.|1998| l. We 
may therefore expect to see measurable differences between 
the dwarf populations in outer regions of clusters (where field 
galaxies are accreted) and the denser central regions, where 
destructive environmental factors also come into play (see 
Faint populations attached to infalling galaxy groups 



7.61 



may also increase dwarf galaxy densities in clusters, although 
they would be more susceptible to tidal stripping than brighter 
group galaxies and may therefore be expected to be located 
preferentially away from the core of the cluster (Biviano et al. 
|1996| l. It is possible to directly search for such environmental 
differences by observing different regions within a cluster and 
constructing region-specific LFs in a consistent fashion. 

At a redshift of z 0.023, Coma ( Abell 1656) is o ne of the 
most extensively studied clusters (see Biviano|1998 for a his- 
torical review), and its richness and proximity make it an ideal 
target for the study of the faint-end LF. It is approximately 
spherically symmetric, and has two central dominant elliptical 
galaxies (NGC 4874 & NGC 4889). Although once thought to 
be fully relaxed, several sub-structures have been found at op- 



tical and X-ray wavelengths (e.g. |Mellier et al. 1988| White et 
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e central 



1993 Neumann et al. 2003 ). 

axy NGC 4874, another on a sub-group of galaxies around 
the giant elliptical NGC 4839, located 40' south-west of the 
center. Recent X-ray observations with XMM-Newton have 
shown the NGC 4839 sub-group to be falling into the cluster 
core ( [Neumann et al.|2001| l. 

An accurate determination of the faint-end slope of the LF 
of the Coma cluster is an essential step towards constraining 
the mass function of galaxies in the local Universe. Numer- 
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ous studies, mainly at optical wavelengths, have measured 
steep faint-end slopes (a ^ —1.7 to —2.3, e .g. Trentham 
[1998) [Milne et al.||2007) [Adami et al.[|2007] l, while others 
see a fai rly flat distribution (e.g. a ~ —1.2, Mobasher et 
[al.|2003| l. These discrepancies are likely, however, to be due 
to both the different photometric limits of the surveys and the 
different methods of subtracting the background galaxy com 
ponent (see 



7.4 1 



In this paper we aim to improve the determination of the 
dwarf galaxy population of Coma by measuring its near-IR 
(3.6/im) LF in two regions of different densities with the In- 
frared Array Camera (IRAC, [Fazio et al.|2004a| on board the 
Spitzer Space Telescope. Observations in the near-IR provide 
us with a different perspective on the cluster LF to that at 
shorter wavelengths, as we observe emission from old stellar 
populations as opposed to the younger populations detected 
in the optical. The near-IR is therefore less sensitive to star- 
formation in galaxies and provides a more direct measure of 
the underlying galaxy mass function of the cluster. Also, near- 
IR measurements are relatively free of biases due to dust ob- 
scuration, which means we have the potential to detect low- 
luminosity dusty galaxies that could be missed in optical sur- 
veys. In the past, the small field-of-view (FOV) of ground- 
based near-IR detectors has made extensive surveys of nearby 
clusters difficult. However, IRAC efficiently mosaics large ar- 
eas of the sky, and has both sufficient spatial resolution (~ 2") 
and sensitivity to perform such a survey. In the shortest wave- 
length IRAC band (3.6/im), the IR zodiacal light is at a min- 
imum of ~23.5 mag arcsec" dKelsaU et al. 1998f, a factor 



of 10,000 low er than ground-based JC-band imaging (^^13.5 



mag arcsec ,|Wainscoat & Cowie 1992). The sensitivity and 



spatial resolution of the IRAC 3.6/Ltm observations, together 
with a redshift-based statistical background source subtrac- 
tion, foreground star subtraction based on Galaxy star counts, 
and conservative optical color filters to guard against back- 
ground galaxies, will allow us to make a better determination 
of the size of the dwarf galaxy population of Coma. A well 
established 3.6/im LF at z can then be compared with 
rest-frame 3.6/im LFs at higher redshifts derived from other 
IRAC channels. For example, 8/im IRAC observations sam- 
ple 3.6/im emission at z ^ 1, which means we could study 
the equivalent near-IR LF (i.e. the mass function of galaxies) 
at this redshift to see how massive clusters evolve. 

In §|2]we describe the Spitzer IRAC observations and data 
reduction procedures. In § [3] we describe our source de- 
tection and photometry, and methods used to separate stars 
from galaxies. In § |4[ we outline the optical photomet- 
ric/spectroscopic surveys we use to identify Coma member 
galaxies, and in § [5] we determine what fraction of IRAC de- 
tected galaxies are members of Coma. In § [6[ we construct 
IRAC LFs for both the core and off-center fields. The results 
are discussed in § [7] and compared with those found at other 
wavelengths. Conclusions are given in § [S] Throughout this 
paper we use a distance modulus (DM) of m — AI — 35.0 
mag, corresponding to a distance to Coma of 100 Mpc for 
iJo=70kms-i Mpc-i. 

2. OBSERVATIONS AND DATA REDUCTION 

The Coma cluster is approximately spherically symmetric, 
with a virial radius of 100' (2.9 Mpc) and total mas s, esti- 



mated from dynamical modeling, of 1.4x lO^^/iyQ^M© dLokas 



& Mamon''2603). The aim of this study is to characterize 




the near-IR dwarf galaxy population of Coma, and to inves- 
tigate whether the LF changes, due to e.g. accretion of field 



1 =1.75 Mpc 



Fig. 1 . — DSSI image of the Coma cluster showing the location of the two 
Spitzer IRAC fields (solid boxes). The optical spectroscopic survey fields 
of Mobasher et al. (2001 1 are shown with dashed boxes, and the circle de - 
notes the virial radius of the cluster (100' or 2.9 Mpc, iLokas & Mamonl2003l. 
The three giant elliptical galaxies are also marked; NGC 4889, NGC 4874 & 
NGC 4839. There are 420 optically identified member galaxies in these two 
fields. 



dwarf galaxies/environmental factors, between regions of dif- 
ferent cluster-centric radius and density. We select two re- 
gions of Coma to overlap with the spectroscopic survey per- 
formed by [Mobasher et al.| ( [200 l| l; these are shown in Fig- 
ure[T] The optical fields are each ~ 32J5x50.'8; one is located 
at the core (Coma 1, centered at 12''59"'23!7, H-28°01'12'.'5), 
and the second south-west of the core (Coma 3, centered 
at 12''57"07!5, H-27°ll'13'.'0). These two fields were cho- 
sen firstly because of their large density contrasts (the local 
density of Coma 3 is somewhere between Coma 1 and the 
field), and secondly because the Coma 3 field covers the X- 
ray secondary peak (the giant elliptical galaxy NGC 4839). 
We refer to Coma 1 and Coma 3 as the 'core' and 'off- 
center' fields respectively. The center of the Coma 3 field 
is located ~57' (1.7 Mpc) from the established central galaxy 
NGC 4874 ( [Kent & Gunn[[T982l l, i.e. ~0.6 times the virial 
radius. 

2.1. Spitzer IRAC Observations 

The Coma 3 off-center field was observed by Spitzer on 
June 13, 2005 (Program ID 3521). The observation consists 
of two pointings designed to cover as much of the Coma 
3 field as possible. AOR 11067904 (Top) is centered on 
12''57'"25!0, H-27°22'13'.'0, while AOR 11067648 (Bottom) 
is centered on 12''57™25!0, H-26°58'13'.'0, the two fields over- 
lapping in the middle. Each pointing is made up of an 8x9 
grid of 72 s per pixel exposures, dithered over 6 positions 
with a medium scale cycling pattern to fully sample the point 
spread function (PSF). Since there are a few moderately bright 
sources in the field, the high dynamic range full array mode 
was used, which includes some very short exposures (0.6 s) 
to aid photometry and PSF fitting of saturated sources. Note 
that these sh ort ex posures are not used when combining the 
images (see § 2.2 1, as this would result in higher background 
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Fig. 2. — The 3.6/jm Spitzer IRAC mosaic of the core (Coma 1, 0.733 
deg^) and off-center (Coma 3, 0.555 deg^) regions of the Coma cluster. 
The fields of the optical photometric/spectroscopic survey of [Mobasher et| 
|al.H200l) are marked with rectangles, each 32.5x50.8'. The circles denote 
the 410 spectroscopically confirmed members of the cluster detected above 
3(T in the IRAC data. 



noise and decreased sensitivity. 

To study the Coma 1 core field, we use archival data taken 
on the core of the Coma cluster. These data, gathered by 
the IRAC instrument team in July 2004 (Program ID 25, 
RI. Fazio), are centered on 12''59'"48?7, +27°58'50'.'0 and 
are slightly deeper than the Coma 3 off-center data (90 s per 
pixel). This observation is comprised of two 12x6 grids 
(AORs 3859712 & 3859968) taken in standard full aiTay 
mode, dithered over 3 positions with a small cycling pattern. 

2.2. IRAC Data Reduction 

The Basic Calibrated Data (BCD) were supplied by the 
Spitzer Science Center (SSC), consisting of individual cali- 
brated frames. The pipeline processing steps performed to 
produce these data were dark subtraction, flat fielding and flux 
calibration (see the IRAC Data Handbook for further details). 
We performed post-BCD processing on the GO and archival 
data in the following manner The individual images were 
mosaiced together using the SSC's data processing package, 
MOPEX, specifying a standard pixel size of 1.22". In order to 
exactly align the images of the 4 IRAC bands. Fiducial Im- 
age Frames (FIFs) were created for each observation using all 
BCD images in that observation and applied to the mosaicing 
of each individual band image. A full mosaic of the two fields 
is shown in Figure |2] 

The images are flux calibrated in the SSC pipeline in sur- 
face brightness units (MJy sr^^). For the purposes of deter- 
mining individual source fluxes (see § |3.1| l, the images were 
converted to flux density units by multiplying by a conversion 



factor of IMJy sr^^ = 34.98 /xJy per pixel. This is a combi- 
nation of the conversion given in the Spitzer Space Telescope 
Observers Manual (IMJy sr^^ = 23.50443 ^Jy arcsec"^), and 
the pixel solid angle (1.22 arcsec^, since the flux density of an 
object is the integral of the surface brightness over its solid 
angle). 

3. SOURCE DETECTION, PHOTOMETRY & STAR-GALAXY 
SEPARATION 

3.1. IRAC Source Detection & Photometry 

Source detection was carried out with the SExtractor (SE) 
algorithm ( Berlin & Arnouts|1996| l. The SE detection thresh- 
old was set to 3(t, with the SEEING_FWHM parameter fixed 
to 1.7" to match the IRAC 3.6/ im PSF full -width-at-half- 
maximum (FWHM, see table 3, Fazio et al.' 2004a). Only 
sources within the areas observed in all four IRAC bands were 
included, to be able to calculate IRAC colors for every source. 
We used the 'double-image' mode in SE, detecting objects in 
the image with the highest signal-to-noise (3.6/im). This pro- 
vides photometry on 4 IRAC bands, and allows a consistent 
measure of IRAC colors. The total areas analyzed in each 
field were 0.733 deg^ (Coma 1) and 0.555 deg^ (Coma 3). 

We measured 'total' fluxes using flexible elliptical apertures 
(i.e. FLUXJ\UTO in SE), whose elongation and position an- 
gle are defined by the object's light distribution, and therefore 



proportional to the size of each galaxy in the image ( [Bertin & 



Arnouts|1996| ). The Spitzer SWIRE team has found fliat fliis 



is the best approach for measuring the flux from small (a few 
arcsecond size) extended sources (see section 5.8 in the IRAC 
Data Handbook v.3). The AUTO elliptical apertures were set 
to 2.5 times the characteristic 'Kron' rad ius of each source, 
so that they contain > 95% of source flux ( Kron|198(){ l. Since 
the Coma field is dense, the local background was subtracted 
using a 5 pixel radius annulus around each source to mitigate 
contamination from nearby sources. The measured flux den- 
sities were then converted to apparent magnitudes using the 
flux density zero points (relative to Vega) quoted in [Reach] 
eTall ( |2005) ; these are 280.9 Jy (3.6Mm), 179.7 Jy (4.5/im), 
1 15.0 Jy (5.8/im) and 64.13 Jy (8/im). Finally, the images and 
source detections were visually inspected. Spurious sources 
(e.g. those in the trails of very bright stars) were removed. 

The IRAC catalog, which will be published separately, con- 
tains four-band data for a total of 29,208 sources, 17,872 in 
Coma 1 and 11,336 in Coma 3. Figure [3] shows the flux den- 
sity distribution of the sources detected in both fields. These 
turn over due to the onset of incompleteness at ~ 15/iJy 
(m^l8.2) in the Coma 1 core field and ^ 21/iJy (m'^17.8) 
in the Coma 3 off-center field at 3.6/im (shown with vertical 
lines). However, to guard against any minor effects of incom- 
pleteness at brighter magnitudes, we adopt more conservative 
completeness limits of ma g^rjx = 17.0 and 16.5 (correspond- 
ing to M3 6;^m = -18 and -18.5) for Coma 1 and 3 respectively 
for our LF calculations (see §[6]). 

To give an insight into the different galaxy populations 
in Coma, Figure [4] shows an IRAC color-color diagram 
(5.8/3.6/im vs 8.0/4.5/im) for the spectroscopically confirmed 
bright (ms.g^m < 14.5) and faint (14.5 < ma.e^m < 18) 
members in the two fields (see § |5] for member identifica- 
tion). This combination of colors is a good indicator of 
star-formation activity in galaxies (i.e. ratios of PAH fea- 
tures to stellar light). The four areas denote the expected 
colors for PAH-dominated (star-forming), stellar-dominated 
(passive) and continuum-dominated (active; AGN) galaxies in 
the four IRAC bands at zero redshift, as predicted by Sajina et 
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Fig. 3. — The 3.6^m flux density distribution of sources detected in the 
IRAC data. The open histogram represents the distribution of the 17,872 
sources detected in the Coma 1 core field, and the shaded histogram shows the 
same for the 1 1336 sources in the Coma 3 off-center field. The solid vertical 
and dashed lines show where major incompleteness starts to affect the Coma 
1 data (15 /iJy, m-iaij.m. ~18.2) and Coma 3 data (21 ^Jy, ms.Sfim "^17. 8) 
respectively. To guard against any minor incompleteness effects brighter than 
this, we adopt more conservative limits of ms.sjim = 17.0 for Coma 1 and 
ms.e^m = 16.5 for Coma 3 for the LF calculations. 



al. (2005 1. This shows that the majority of the Coma members 
have colors consistent with those of passive stellar galaxies, 
as expected for such a massive cluster Another population 
have colors consistent with AGN, which can be confirmed 
by follow-up X-ray observations with Chandra and XMM- 
Newton. Note that the bright members are skewed towards 
the passive galaxy colors, while the faint members randomly 
populate the color-color space. 

3.2. Star-Galaxy Separation 

According to the IRAC source number counts of Fazio et al.' 
( [2004 b ), a substantial fraction of sources in the IRAC bands, 
particularly at bright magnitudes, are foreground Galactic 
stars. In order to reliably measure the number of galaxies 
in the data, it is important to filter these out. At the brightest 
magnitudes in our catalog (m3 g^m < 14.5), we can use the 
spatial extent of galaxies at 3.6//m to morphologically sep- 
ar ate them from stars. He re we use the me t hod de scribed 
in [Fazio et al.] P004b| l and [Eisenhardt eTal^ ( [2004l i, which 
uses the difference between small 3" diameter circular aper- 
ture and 'auto' aperture flux measurements in SB as a 'con- 
centration parame ter' . In this scheme , which is reliable down 
to m3 6;^m ~14.5 (Fazio et al. 2004b I, point sources will have 
approximately the same magnitude using both 'auto' and cir- 
cular (corrected) apertures, whereas galaxies will show an ex- 
cess in the corrected circular aperture because their emission 
does not match the IRAC PSF profile. 

Figure |5] illustrates the concentration parameter, showing 
the boundaries defined by jEisenhardt et aL[p004) to separate 
stars from galaxies. Sources with values of the SE parameter 
CLASS_STAR > 0.95 ai'e also shown (where values close to 
1 represent point sources), demonstrating that these sources 
do indeed fall inside the designated ranges. An aperture cor- 
rection of -0.69 mag was derived from the data to correct the 
3" aperture magnitudes out to the IRAC total magnitude cal- 
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Fig. 4. — IRAC color-color diagram of the confirmed Coma members 
detected in all 4 bands in both fields, defined in terms of the logarithms of 
flux densities in the 5.8/3.6/.tm and 8.0/4.5/jm bands. The dashed lines mark 
boundaries between different galaxy types as defined by Sajina et al. 1 2005J. 
Green diamonds and red squares denote the blight (ma.opm < 14.5) spec- 
troscopic members in Coma 1 and 3 respectively, while the fainter members 
are shown with crosses. The majority of the bright Coma members have col- 
ors consistent with passive stellar galaxies. 



ibration diameter of 24.4" (20 pixels). This is defined, in 
terms of flux, as the average ratio (1.90) of the 24.4" to 3" 
aperture measurements from SE on 56 unsaturated stars (9 
< m3 < 12) using a 5 pixel-wide background annulus 
(this magnitude range excludes the few brightest stars that 
are obviously saturated). Note that this small aperture size 
does not provide a reliable measurement of the total fluxes of 
galaxies, but it is useful for the purpose of star-galaxy sep- 
aration. All sources satisfying the criterion -0.25 < ma g^m 
(3") - m3.6/jm (auto) < 0.2 and brighter than 14.5 mag were 
flagged as possible stars. These were all inspected visually, 
and 7 were re-classified as galaxies based on their IRAC mor- 
phology. As an additional check, the positions of the sources 
classified as stars were cross-matched with the optical photo- 
metric catalog oflKomiyama et al. ( 2002) (which has also had 
the stars removed by a similar method) using a 3a matching 
radius of 4.24" (see §[5]l, and only 4 sources were reclassified 
as galaxies based on these matches. The 4 brightest stars in 
these fields were also classified as galaxies by their concentra- 
tion parameter due to the bright wings of their PSFs, and these 
were subsequently re-classified as stars after visual inspection 
of the data. In total, 476 bright sources (ms g^m < 14.5) were 
classified as stars over the two fields, and removed from the 
LF calculations. We also removed 2100 sources with ma g^jm 
(3") - ma. (auto) <-0.25 fainter than 14.5 mag as these 
are almost certainly cosmic rays ( [Eisenhardt e t al. 2004). 

At fainter magnitudes (ma.e^m ^ 14.5), it is difficult to dis- 
tinguish stars from galaxies morphologically, so we must sub- 
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tract the stellar contribution using statistical methods. Here 
we use the "DIRBE (Diffuse Infrared Background Experi- 
ment) Faint Source Model" for Galactic stars, which is an im- 



(auto) 



Fig. 5. — The stellar concentration parameter, as defined by [Fazio et al.| 
j20Q4b) andlEisenhardt et al. (20041 (see text for details). AH sources de- 
tected in both Coma fields are marked with small points, and those with the 
SE CLASS-STAR parameter > 0.95 are over-plotted with stars. The hori- 
zontal lines denote the boundaries of the stellar classification, and the vertical 
line shows the magnitude down to which this is reliable (ms gj^m = 14.5). 
Sources below the lower bound of the stellar classification and fainter than 
iTi3.6nm = 14.5 (lower right-hand comer) are classed as cosmic rays (CRs). 
All sources brighter than ms.Sfim = 14.5 classed as stars, plus fainter 
sources classed as CRs, are excluded from the LF calculations (see text for 
further details). Fainter than this, stars are subtracted from the LFs using the 
DIRBE Faint Source Model (see Figure[6|. 
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Fig. 6. — Differential number counts of all sources detected at 3.6/.tm in 
the Coma 1 core field {green diamonds) the and Coma 3 off-center field (red 
squares). The predi cted number coun ts of Galactic stars from the DIRBE 
Faint Source Model jArendt et al.|I998j for the line-of-sight towards Coma 
(Galactic coordinates: L = 58.U8, b = 87.96) are shown as a dashed line. 
At bright magnitudes {^z.fi^m < 14.5), stars are subtracted from the LF 
calculations using the stellar concentration parameter (Figure[5l. Fainter than 
this, the stellar contribution is statistically subtracted using the Faint Source 
Model. 



plementation of the Wainscoat et al. ( 1992) statistical model 
(further details can be found in Arendt et al.| |1998). Fig- 
ure |6] shows the differential number counts of the two Coma 
datasets, together with the predicted model counts for Galac- 
tic stars in the DIRBE 3.5/im band for the line-of-sight to- 
wards Coma (Galactic coordinates: I — 58.08, b — 87.96, 
R. Arendt, private communication). The model demonstrates 
the dominance of stars at bright magnitudes (8 < ms.gpm < 
14). However, the counts in the central Coma 1 field are well 
above the star counts from as bright as ma g^,„ 10 due to 
the high number of bright galaxies located there, whereas the 
counts in the off-center field do not diverge significantly from 
the model until m3 g^m ~ 13. Fainter than m3 fi^„i ^ 14.5, 
there is a much smaller stellar contribution to the total num- 
ber counts (^10-15%), and we subtract the model star counts 
from our data in this range for our LF calculations in § |6] 

3.3. Comparison with field galaxy counts 

In Figure [7] (a) we plot the Coma galaxy number counts 
(excluding all contribution s from stars) together with the field 
galaxy number counts of |Fazio et al. (2004bj). These con- 
sist of observations from 3 fields: QSO 1700 (5' x 10'), Ex- 
tended Groth Strip (EGS) (0?17 x 2°) and Bootes (3° x 3°). 
For magnitudes brighter than 15 mag, we use the wide-area 
Bootes field counts; fainter than this we use the average of 
the counts from the d eeper QSO 1700 an d EGS survey fields 
(taken from table 1 in Fazio et al.|2004b| l. This figure shows 
that there is a clear excess of galaxies in both Coma fields at 
bright magnitudes (ma g^m < 14.5) as expected , but it also 
indicates an excess at fainter magnitudes. To illustrate this 
more clearly. Figure |7](b) shows the ratio of galaxy counts in 
the two Coma fields to the field counts in the range 14.5 < 
m3.6/im < 17. Even taking possible cosmic field-to-field vari- 
ati ons into account (est imated as 16-19% over 1100 arcmin^ 
by Dahlen et al. [2005 similar to the EGS deep survey area 
of 1224 arcmin^), the excess is still above the highest back- 
ground count estimate. Note, however, that since we know 
that there is significant background structure behind Coma, 
we do not use statistical background source subtraction for 
the LFs in this study. Instead, we use spectroscopic member- 
ship fractions together with optical color selection to subtract 
the background galaxy component (see §|5]& 



6.1 



4. OPTICAL DATA 

To determine the contribution to the overall LF from back- 
ground galaxies, we use redshift information from two optical 
cat alogs. The first is the spectroscopic survey of Mo basher et| 
aL]( [2fe) , based on the William Herschel Telescope (WHT) 
wide-field photometric B- and i?-band CCD survey of five 
fields in the Coma cluster by |Komiyama et al.| ( |2002[ ). In 
brief, two of the photometric survey fields, each 32.'5x50f8, 
were chosen for follow-up spectroscopic observations using 
the WYFFOS multi-fiber spectrograph on the WHT; Coma 
1 (core), centered on 12''59"23!7, H-28°01'12'.'5, and Coma 
3 (off-center), centered on 12''57™07!5, H-27°ll'13'.'0. In 
this survey, redshifts were measured for two classes of galax- 
ies in Coma: the "bright" sample, consisting of galaxies al- 
ready confirmed as spectroscopic members with previously 
measured redshifts, and a "faint" sample of galaxies with 
no known redshifts, which were on average fainter than the 
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Fig. 7. — (a) Differential number counts of the IRAC 3.6^m sources 
classed as galaxies in the Coma 1 core field (green diamonds) and the Coma 
3 off-center field (red squares). Sources classified as stars or cosmic rays 
have been removed using the stellar concentration parameter and the DIRBE 
Faint Source Model. The dashed line shows the field galaxy counts of [Fazio | 
|eFaL](2004bt ; a clear excess of galaxies is evident in both Coma fields with 
respect to the field population at bright magnitudes, (b) T he rati o of galax ies 
detected in the Coma data to the field galaxy counts of Fazio et al. 1 2004b I for 
the faint magnitude range 14.5 < 1x13 g^jm < 17, demonstrating that there is 
also an excess of faint galaxies in Coma with respect to the field. 
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Fig. 8. — (a) The i?.-band magnitude distribution of the optical spectro- 
scopic targets in the combined Mobasher et al. (2 001^ and|Colless & Durin] 
jl99 6) catalogs for the two Coma fields. The total number oTtarget galaxies 
per 0.5 mag bin is shown with a shaded histogram, while the number of IRAC 
detected galaxies matched to the i?.-band sample is overlaid with an open his- 
togram. Almost all (97%) of the /?-band sample is detected at 3.6/im. (b) 
The 3.6/jm magnitude distribution of the spectroscopically confirmed mem- 
bers in the central Coma 1 data (solid line) and the Coma 3 off-center data 
(dashed line). 



bright sample. The faint sample was chosen from the photo- 
metric data using their positions in the color-magnitude plane 
(B — R) versus R with the criterion 17.5 < i? < 20 (designed 
to include the dwarf galaxy population) and 1 < B — R <2 
(to eliminate the majority of background galaxies). 

We supplement this with a wider-field (2.63 deg'^) but shal- 
lower (i? ~ 18) catalog compiled by M. Colless (private 
communication), to ensure the inclusion of previously-studied 
brighter galaxies in the cluster that were not re-measured in 
the Mobasher et al. ( 20011 study. This supplemental catalog 
includes redshift measurements by |Colless & Dunn| ( |1996 1 
using the Hydra multi-fiber spectrograph at KPNO, measure- 
ments from the Two-Degree Field (2dF) survey ofjEdwards£t] 
aL|(|2002||, and supplemental redshifts from the NASA/IPAC 



Extragalactic Database (N ED) and literature ([ Kent & Gunn 
T982 Cald well et al.|1993|[varrHaarlem et al . 1993: Biviano 
|et al._1995) . In total, these combined catalogs provide spec- 
troscopic redshift measurements of 727 galaxies across the 
Coma 1 and Coma 3 fields. 

Note that due to the color selection of the spectroscopic 
sample, we also use deeper photometric data (r ^ 22.5 mag) 
from the recent Sloan Digital Sky Survey (SDSS) data re- 
lease 5 (DR5) to impose a similar optical color cut (SDSS 
g — r < 1.3) on our IR AC s ample when constructing the 
3.6/j,m faint-end LPs (see §|6.1|l. 



5. IRAC CLUSTER MEMBERSHIP & SPECTROSCOPIC 
MEMBERSHIP FRACTIONS 
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TABLE 1 

Cluster membership fractions for Coma 1, Coma 3 and the total 
of the two 



n^3.6/im 


Coma 1 core 


Coma 3 off-center 




Total 


f(m) 


A f(m) 


t(m) 


A f(m)\ 


/(m) 


A L/(m)J 


9.00 


1.00 


0.00 


1.00 


0.00 


1.00 


0.00 


10.25 


1.00 


0.00 


1.00 


0.00 


1.00 


0.00 


10.75 


1.00 


0.00 


1.00 


0.00 


1.00 


0.00 


11.25 


1.00 


0.00 


0.80 


0.18 


0.95 


0.05 


11.75 


1.00 


0.00 


0.91 


0.09 


0.97 


0.03 


12.25 


1.00 


0.00 


1.00 


0.00 


1.00 


0.00 


12.75 


0.81 


0.09 


0.60 


0.15 


0.74 


0.08 


13.25 


0.79 


0.08 


0.75 


0.15 


0.78 


0.07 


13.75 


0.73 


0.07 


0.50 


0.12 


0.65 


0.06 


14.25 


0.49 


0.06 


0.32 


0.11 


0.45 


0.05 


14.75 


0.59 


0.07 


0.37 


0.09 


0.52 


0.05 


15.25 


0.44 


0.07 


0.37 


0.08 


0.41 


0.05 


15.75 


0.40 


0.07 


0.39 


0.09 


0.40 


0.05 


16.25 


0.50 


0.09 


0.42 


0.11 


0.47 


0.07 


16.75 


0.57 


0.09 


0.25 


0.10 


0.44 


0.07 


17.25 


0.27 


0.11 


0.20 


0.18 


0.25 


0.10 


17.75 


0.33 


0.27 











o 



S3 



0.8 



0.6 



0.4 



0.2 



Coma 1 
Coma 3 
Total 




In order to identify cluster membership, we cross-correlated 
the IRAC catalogs with the two optical Coma spectroscopic 
catalogs discussed in §|4] The positions of both the IRAC and 
optical catalogs are good to within < 1", and a 3(t search 
radius of 4.24" was used for the positional matching, cor- 
responding to three times the two positional errors added in 
quadrature. Figure [8] (a) shows the i?-band magnitude dis- 
tribution of the spectroscopic targets in the Coma 1 core and 
Coma 3 off-center fields, plus the numbers of those galax- 
ies detected in the IRAC 3.6/im data. The detection rates 
with IRAC are high; out of a total of 498 spectroscopic target 
galaxies in the Coma 1 core field, 487 (98%) were detected 
by IRAC, while in the Coma 3 off-center field, 221 (97%) out 
of 229 were detected. 

In the spectroscopic catalogs, galaxies with recession ve- 
locities in the range 4000 km s~^ < cz < 10,000 km s^^ are 
considered to be members of Coma, corresponding to the 3cr 
range measured by Colless & Dunn ( 1996). In the Coma 1 
field, 312 of the galaxies are members, 306 (98%) of which 
were detected with IRAC. In Coma 3, 108 galaxies are mem- 
bers, of which 104 (96%) had IRAC matches. The positions of 
these 410 confirmed members are plotted on the IRAC 3.6 /im 
mosaic in Figure |2] and their 3.6/im magnitude distributions 
are shown in Figure [8] (b). IRAC finds almost all of the known 
optical members. Note that there is some weak IR emission 
at the locations of the few undetected members in the 3.6/im 
images (all fainter than R= 16.5), but not significant enough 
to be detected above the 3fT threshold. 

To construct LFs for the Coma galaxies, we must estimate 
the contribution of background galaxies to the IRAC number 
counts. To achieve this, we compute a cluster membership 
fraction for each field. The opt ical spectroscopic cover age of 
Coma is not complete, but the 'Mobasher et al. ( 2003| spec- 
troscopic survey was specifically designed to sample optical 
color magnitude space in a way that determines membership 
within each half-magnitude range down to i? ^ 20. Thus, al- 
though not every galaxy down to i? = 20 has a spectrum, the 
membership of the cluster is constrained to this depth. There- 
fore, we assume that the incomplete spectroscopic sample is 
representative of the whole galaxy population in the Coma 
cluster, so that the fraction of galaxies determined to be cluster 
members in the spectroscopic survey as a function of magni- 
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Fig. 9. — The spectroscopic membership fractions of the IRAC 3.6/im 
detected galaxies as a function of apparent magnitude for the Coma 1 and 
Coma 3 fields, and the total of the two. 



tude can be applied to the complete catalog of IRAC detected 
galaxies. The magnitude limits of the IRAC and optical spec- 
troscopic catalogs are well matched for determining this func- 
tion: the faintest 3.6/im bins correspond to mj^ ^ 20 — 21, 
the limit of the spectroscopic data. 

We define the spectroscopic membership function (f, 
Eqn. [T} as the ratio of the number of spectroscopically con- 
firmed cluster members in a given apparent magnitude bin, 
A'c(m|spec), to the total number of galaxies in the spec- 
troscopic catalog (members and non-members), detected by 
IRAC, in the same magnitude bin, Nt(m|spec): 



/(to I spec) 



Nc{m I spec) 
Nt{m I spec) 



(1) 



This ratio is calculated separately for the Coma 1 core and 
Coma 3 off-center fields (and total of the two). The relative 
error on these functions, assuming poisson statistics, is given 
by: 



df{m\ spec) 
f{m\ spec) 



Nc{m I spec) Nt{m \ spec) 



(2) 



The membership fractions were determined in half- 
magnitude bins for galaxies with ma > 10. At the bright 
end, where there are far fewer galaxies, they were binned be- 
tween 8 < m3.6;^T?i < 10 to ensure that the statistics were good 
enough for constraining the data. These membership fractions 
(and errors) per bin for the two fields and for the total of the 
two are shown in Table [T] and plotted in Figure |9] The Coma 
1 data extends down to ma.gpm = 18, while the Coma 3 is 
limited to m3 gp,„ = 17.5 due to the slightly shallower obser- 
vation, although we impose slightly brighter magnitude limits 
for the LF calculations (see §|6]l. As expected, the fraction of 
member galaxies per magnitude bin is higher in the Coma 1 
field due to the relatively larger number of brighter galaxies at 
the center of the cluster 
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TABLE 2 

Luminosity functions for Coma 1, Coma 3 and the total of the two 





Coma 1 (core) Coma 3 (off-center) 


Total 




9i(M) A[(/i(M)J (/.(M) A[<^(M)J 


0(M) A l^(M)\ 




All IRAC detections (no optical constraint) 



-26.00 0.56 0.50 

-24.75 5.37 1.55 

-24.25 8.95 2.00 

-23.75 8.06 1.90 

-23.25 12.98 2.41 

-22.75 17.46 2.80 

-22.25 11.23 2.34 

-21.75 17.36 3.08 

-21.25 20.58 3.31 

-20.75 28.48 4.30 

-20.25 50.65 6.73 

-19.75 94.84 15.39 

-19.25 169.38 28.51 

-18.75 368.62 67.91 

-18.25 670.16 110.46 

-20.25 45^3 6.16 

-19.75 75.63 12.40 

-19.25 111.73 18.99 

-18.75 185.54 34.48 

-18.25 235.84 39.37 

9 - r 

-20.25 19.78 37l8 

-19.75 34.85 6.03 

-19.25 63.92 11.09 

-18.75 117.05 21.98 

-18.25 160.89 27.10 



0.44 


0.51 


0.51 


0.36 


0.59 


0.59 


3.31 


0.92 


2.36 


1.18 


6.11 


1.25 


2.36 


1.18 


5.55 


1.20 


6.98 


2.05 


10.34 


1.63 


5.32 


1.77 


12.22 


1.76 


3.90 


1.55 


7.94 


1.49 


7.09 


2.29 


12.93 


2.01 


10.63 


3.07 


16.50 


2.32 


8.02 


2.97 


19.91 


2.78 


31.28 


8.27 


44.20 


5.27 


64.06 


14.58 


81.24 


10.77 


140.85 


33.58 


157.23 


21.62 


308.89 


83.56 


345.32 


52.84 



r < 2 (r < 22.5) faint end 



29.75 
51.56 
100.01 
155.44 



7.89 
11.83 
23.99 
42.27 



40.77 
65.01 
106.81 
173.79 



4.90 
8.70 
14.80 
26.79 



g - r < 1.3 (r < 22.5) faint end 



11.16 
21.72 
50.82 
95.75 



3.21 
5.25 
12.42 
26.22 



16.62 
28.98 
58.32 
108.53 



2.30 
4.09 
8.25 
16.87 



Note. — LFs are in units of galaxies 0.5 mag ^ Mpc ^. A distance modulus 
of 35.0 mag (i?o=70km s^^ Mpc^^) is assumed for Coma. 



6. THE LUMINOSITY FUNCTIONS 

In this section we present the IRAC 3.6/xm LF of galaxies 
in the Coma cluster. The LFs for the Coma 1 and Coma 3 
fields are constructed separately, using the membership frac- 
tions found in the last section. These are estimated as: 



(/)(m) 



Npjm) 
A 



f{m I spec) 



(3) 



where Np{m) are the galaxy counts, / is the membership 
fraction and A is the area of the cluster over which the counts 
are measured. To improve statistics on the faint end of the 
IRAC LF, we extend our sample beyond the regions for which 
spectroscopic data are available. The relative errors on the 
LFs are given by: 



d<j){'m) 
<j){m) 



Np{m) Nc{m \ spec) Nt{m \ spec) 



(4) 

At a distance of 100 Mpc, 1°= 1.75 Mpc, and the areas cov- 
ered by the IRAC observations correspond to 2.234 Mpc^ for 
the Coma 1 core field, and 1 .693 Mpc^ for the Coma 3 off- 
center field. "Raw" LFs for both fields, calculated using all 
IRAC-detected galaxies, are listed in Table |2] and plotted in 
Figure [TO] However, in the next section we impose opti- 
cal color selections to the faint sample to create optically- 
constrained faint-end LFs, and we therefore consider the faint 
end of these raw LFs (Ms.g^m > —20.5) to be upper limits 



only. The Coma 1 LF covers the absolute magnitude range of 
—27 < Ms.g^tm < —18 (using a distance modulus of 35.0), 
while for Coma 3 the LF extends down to M3 g^„j —18.5. 
Within these magnitude ranges, these initial calculations im- 
ply a total of 3320 IRAC-selected members in the Coma 1 
core field and 1006 in the Coma 3 off-center field. This in- 
cludes both the known spectroscopic members and additional 
galaxies with no redshift information (see Table [3]l. The LFs 
for both fields are consistent with a Schechter form, with a 
striking upturn at Ma g/jm ~-20.5, presumably caused by the 
onset of the dwarf galaxy population. 

6.1. Optical constraints 

As mentioned in § |4j the spectroscopic sample used to de- 
fine the membership fraction is selected based on optical color 
criteria. Therefore, we have used optical data from the SDSS 
DR5, which covers both Coma fields, to create two subsets 
of IRAC-detected sources filtered by their optical properties 
(g — r color and r-band magnitude) that are more closely 
matched to the optically-selected population on which spec- 
troscopy was performed. We seek to strike a balance between 
guarding against this optical selection bias and losing valuable 
information afforded by our IR selection. We first require that 
the IRAC sources have optical counterparts above the SDSS 
photometric limit of r ~ 22.5 (i.e., we do not include any op- 
tically undetected sources in our LFs). Inspection of the g — r 
versus r plot in Figure 1 1 shows that the optical spectroscopy 



of Mobasher et al. (2001) samples this color-magnitude space 
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Fig. 10. — The 3.6/im LFs for the Coma 1 core field (top left, green), the Coma 3 off-center field (top right, red) and the total of the two (lower left, black). The 
best-fit Schechter functions, listed in Table[3] are shown at the bright end of each LF (ms.g^tm < 14.5). At the faint end (ms g^^m > 14.5), we show the 'raw' 
LFs, which take into account all IRAC-detected galaxies, and which we consider to be upper limits to the faint-end normalization. We also show the two g — r 
color-selected faint-end LFs for the IRAC sources with SDSS optical counterparts down to r = 22.5, with the g~r < 1.3 selection being the most conservative. 
The faint-end power-law fits, also listed in Table[3] are shown with straight lines. Note that although the overall number of galaxies decreases with these optical 
selections, the faint-end upturn is still present. In the lower right panel, we also show a comparison between different methods of baclcground subtraction for the 
total LF (membership fraction and field subtraction); the two are consistent in almost all bins, and the faint-end upturn is present in the field-subtracted LF. 



to g — r ~ 2. Note that there is a known background clus- 
ter behind Coma at z 0.5 (Adami et al. 2000), whose red 
sequence is clearly visible in Figure 11 around g — r ^ 1.4. 
Using the color-magnitude distribution of galaxies in the local 
Universe found in |Hogg et al. (j2004), we estimate ^ 97% of 
galaxies to have g — r < 1.2. With all these considerations in 
mind, we find that g — r < 1.3 is the best compromise color, 
which we term "conservative", in that it is the most strict color 
filter that does not over-subtract the galaxies but also guards 
against contamination by the background cluster We have 
also included an "optimistic" g — r < 2 sample, correspond- 
ing to the co lor space sampled spectroscopically by |Mobasher| 
|etaL] ( |200T] i. 

Overall, ~60% of IRAC sources above the adopted 3.6/im 
LF completeness limits have optical counterparts down to 
r = 22.5. When the color criteria are applied, there are 4,894 
3.6/xm galaxies in the optimistic sample, and 3,051 galaxies 
in the conservative sample, after subtraction of the faint-end 



stellar contribution (see Table [3]l. We choose to only apply 
the color cuts to the faint end of the LFs (ms g^^m > 14.5), 
where the largest uncertainty in the cluster membership frac- 
tion exists (due to poorer spectroscopic sampling statistics). 
When the spectroscopic membership fractions are applied to 
the g — r < 2 filtered galaxy counts, a significantly reduced 
number of members is implied: 1757 in Coma 1 (—27 < 



M: 



3.6/j.m 



< -18) and 653 in Coma 3 (-27 < M, 



< 



-18.5). For g — r < 1.3, the member totals are further re- 
duced to 1181 in Coma 1 and 387 in Coma 3. The two sets of 
optically filtered LFs are plotted for both fields in Figure 10 
together with the optically unconstrained "raw" IRAC LFs to 
demonstrate the plausible ranges of slope and normalization. 
Note that although the numbers of member galaxies are re- 
duced by the color-cuts, the rise at the faint end is still signif- 
icant in both fields. 

6.2. Parametric Fits 
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TABLE 3 

SCHECHTER AND POWER-LAW FITS TO THE COM A 1, COMA 3 AND TOTAL LFS 



Optical color selection ^3 6m"i Schechter slope (ai) Faint-end slope (02) ^gal -'V^e^je^ (inferred) 

(1) (2) (3) (4) (5) (6) 



Coma 1 (core): -27 OVIg.e^m < -18 



None 


-25.17±0.10 -1.18±0.06 


-2.40±0.13 


6280 


3320 


g-r <2(r < 22.5) 




-1.93±0.13 


3326 


1757 


g-r < 1.3 (r < 22.5) 




-2.18±0.13 


2156 


1181 




Coma 3 (off-center): -27 <M3 


6)im < —18.5 






None 


-24.69±0.15 -1.30±0.08 


-2.68±0.30 


2425 


1006 


g-r <2(r < 22.5) 




-2.23±0.30 


1568 


653 


g - r < 1.3 (r < 22.5) 




-2.60±0.33 


895 


387 




Total: -27 <M3.6^m < 


-18.5 






None 


-25.28±0.20 -1.25±0.05 


-2.48±0.15 


6086 


2846 


g-r <2{r < 22.5) 




-2.05±0.15 


3973 


1898 


g - r < 1.3 (r < 22.5) 




-2.38±0.15 


2423 


1215 



Note. — (1) SDSS optical color selection of the faint (M3.6;jm > —20.5) IRAC 3.6fim galaxies. "None" 
refers to IRAC LFs with no optical constraint. (2) Bright magnitude turnover in the Schechter function. (3) 
Schechter slope, Ofi, of bright end of the LF (M3. g^im < —20.5). (4) Faint-end slope (Ms.efim > —20.5), 
where 02 = — (fc/0.4 + 1) and k is the logarithmic power-law slope. (5) Total numbers of galaxies detected 
at 3.6/im over given magnitude ranges, excluding all stars (morphological/statistical). (6) Total numbers of 
member galaxies inferred by the spectroscopic membership fractions and optical color-cuts over the given 
magnitude ranges. 
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Fig. 1 1. — The g-r versus r color-magnitude relation for the IRAC 
3.6/^m detected galaxies matched to optical sources with r < 22.5 in the 
SDSS DR5 catalog. The Coma member galaxies form a strong red sequence 
at0.5<g — r<l, and also note the red sequence of the background cluster 
at g — r ~ 1.4. The dashed horizontal lines denote the two optical color cuts 
applied to the faint end of the LF. The g — r<2 cut is used to coiTect for any 
optical bias in the spectroscopic completeness function, and corresponds to 
the color space sampled by the Mobasher et al. 1 2001 1 spectroscopic survey. 
The g-r < 1.3 cut is our most conservative, and is chosen to subtract the 
contribution of the background cluster and to be consistent wi th the SDSS 
color- magnitude distribution of galaxies in the local Universe (Hogg et al.| 
[20041 . 



The widely used parametric form to describe the shape of 
the galaxy LF is the Schechter ftinction ( Schechter|1976[ ): 



$(M) = 0*X"+ie^^ (5) 

where X = iq-O-^m-m') ^ characteristic mag- 

nitude where the distribution turns over at the bright end due 
to the sharp decline in the numbers of galaxies at these mag- 
nitudes, a is the logarithmic faint-end slope, and 0* is the 
normalization. 

A lot of studies of the Coma LF fit a single Schechter func- 
tion to the data. However, the sharp rise at the faint end of the 
3.6/im LF in both fields makes it obvious that a single func- 
tion will not be adequate in this case. Instead, we choose to fit 
a Schechter function to the bright end (Ma e^m < —20.5) to 
model the M* and a parameters, and an exponential power- 
law function with the following form to model the dwarf 
galaxy population at the faint end (Ma.g/im > —20.5); 



kM 



(6) 



where k is the power-law slope, which can be transformed 
into the equivalent o f the Schechter function slope via a = 
-(fc/0.4 + 1) (e.g. [Sandage et"aL]|l985f To distinguish 
between the slopes of the two components, we refer the 
Schechter slope as ai and the faint-end power-law slope as 

The results of the fits are shown in Table |3] for the two 



Coma fields (and their total), and are plotted in Figure 10 
At the bright end, the best-fitting Schechter slopes range be- 
tween fairly flat values of ai = -1.18±0.06 (Coma 1) and - 
1.30±0.08 (Coma 3), with M3 values ranging between 
-25.17±0.10 and -24.69±0.15 respectively. At the faint end, 
the raw IRAC LFs are very steep with slopes 0:2 = -2.40±0.13 
(Coma 1) and -2.68±0.30 (Coma 3). With the most conser- 
vative optical color-cut imposed (g — r < 1.3), these are re- 
duced to a2 = -2.18±0.13 in Coma 1 and -2.60±0.33 in Coma 
3. Note that there is a difference in all these parameters be- 
tween the Coma 1 core and Coma 3 off-c enter fields, which is 
indicative of environmental effects (see 



7.61. 



Spitzer IRAC LF of the Coma Cluster 



11 



7. DISCUSSION 

This study provides, for the first time, a measure of the 
3.6/xm LF for a cluster at z ^ 0, and indicates that IRAC 
observations may provide a very efficient way to search for 
dwarf galaxies in clusters. It covers a similar absolute mag- 
nitude range as IR LFs for galaxies in more distant clusters, 
allowing study of the evolution of the mass function in the 
cluster environment. 

The shape of the IRAC 3.6fim LF at the bright end 
(Ms.g/jTjj < —20.5) is consistent with a Schechter form. The 
characteristic magnitude (M3 g^„) is brighter (by ^ 0.5 mag) 
for galaxies at the denser core region compared to the off- 
center population. Assuming that the 3.6/im flux measures 
the underlying mass of galaxies, we conclude that the con- 
centration of massive galaxies at the core of the Coma cluster 
is presumably caused by mergers of smaller galaxies as they 
fall in the gravitational field of Coma. This is supported by 
the increase in space density of the less massive galaxies (i.e. 
the Schechter faint-end slope, ai) towards the outskirts of the 
cluster 

The important result from this study is depicted in Fig- 
ures [TO] and 13 where we explore the very faint tail of the 
3.6fim LF (a2)- The steep upturn at Ma e^m ^ —20.5 is 
very striking, and is present in bot h the core and off-center 
regions of the cluster. As shown in § |6.1| a significant fraction 
(^--^60%) of IRAC-detected galaxies have optical counterparts 
in the SDSS data down to r ~ 22.5, which indicates that 
this sub-sample are not spurious data artifacts. Using spec- 
troscopic identification and conservative optical color cuts to 
remove contamination by background objects, the upturn at 
the faint end remains. Furthermore, the background cluster at 
z 0.5 is accounted for and removed by the g — r < 1.3 
color selection. Also, an independent technique for back- 
ground subtractio n, us ing counts from field samples, gives the 
same trend (see § |7.1| l. 

The results from the conservative g — r < 1.3 color selected 
IRAC sample infer that there are at least ^ 1600 member 
galaxies over the two fields (see Table 3}. A simple extrapola- 
tion of these results also implies that there are at least ^ 5300 
member galaxies over the whole area of Coma out to the ex- 
treme radius of the Coma 3 field (~ 87'), assuming the same 
number density of galaxies as the Coma 3 field across the un- 
observed area. The inferred stellar masses of the galaxies in 
t he up turn portion of the LF is ~ 5 x 10^ — 3 x 10^ Mq (see 
§ |7.4| i. This is comparabl e to the stellar m ass of, for instance, 
the SMC (~ 3 X 10*^ Mq 
than the LMC (~ 2 x 10^ 



Lero y et al. 2007 J, but less m assive 
^ 



Schommer et al. 



1992: 



However, there are a number of selection effects which 
might lead to this upturn feature, with the most likely being 
contribution by background objects or cosmic variance. Bear- 
ing in mind these caveats, and assuming the observed upturn 
to be real, we now investigate the implications of this result 
by making a detailed comparison with other multi-wavelength 
studies in the literature. Studies of the Coma LF over the past 
~25 years are numerous, particularly at optical wavelengths, 
and we have compiled measurements of the LF slope from the 
literature in Table 4] Also listed are the areas covered by each 
study and the method of background source subtraction used. 
Where single functions were fitted to the data, the value of a 
reflects the slope over the entire magnitude range above M* in 
that particular band. In cases where only the faint data were 
fitted, or where two component models were used (either a 
second Schechter function or a power-law), the quoted value 



of a is for the faint component only. 

7.1. Membership & Background Subtraction 

It is important to consider the method of background source 
subtraction when comparing results, as different methods can 
yield different values of a. Traditionally, three main types 
of background subtraction have been used: (1) statistical sub- 
traction of control field galaxy counts, (2) color information 
(B — R) to exclude galaxies redder than the average for the 
cluster and (3) redshift/membership fraction, the method used 
in this study. Since the advent of wide-field CCD detectors, 
many studies have used nearby offset control fields for back- 
ground source subtraction. However unless the control fields 
are very large (at least the size of a cosmic void), cosmic 
field-to-field variations in background galaxy counts have the 
potential to introduce large uncertainties. This method is re- 
quired for background subtraction at very faint magnitudes 
(R > 22) due to current limitations of multi-object spec- 
troscopy. However, it may not be appropriate for studying 
a rich cluster with a large angular size like Coma, particularly 
when one might expect the entire environment (i.e. the region 
of space behind Coma) to perhaps show over-enhancements 
in galaxy counts due to background groups, filaments etc. In- 
deed, Adami et al.| (j2000) detect three galaxy groups (z ^ 
0.1-0.3) and a more distant cluster at z ^ 0.5 along the line- 
of-sight to the Coma cluster, which would contaminate Coma 
galaxy counts. Using simulations, Valo tto et al.| (|2001) have 
also demonstrated that such projection effects, resulting from 
large-scale structure behind clusters, can produce artificially 
steep faint-end slopes (a < —1.5). 

In this study, we have spectroscopic information for the en- 
tire 3.6/im magnitude range of the LF, which gives a more 
reliable estimate of the faint-end slope in our data. For com- 
parison, however, we have constructed a field-subtracted LF 
using the field counts of (Fazio et al. ( |2004b| l, which is shown 
in Figure 10 (lower right). Although not well constrained, it 
is consistent in almost all bins with the LF computed using the 
spectroscopic membership fraction within the combined pois- 
sonian errors. The faint end of the field-subtracted LF will be 
very uncertain due to the smaller size of the EGS control field 
(1224 arcmin^), but importantly we demonstrate that the steep 
upturn is present in the data regardless of the background sub- 
traction method. 

7.2. Comparison to other Coma IR LFs 

The majority of Coma L F stu dies have been carried out at 
optical wavelengths (see § 7.4 1. However, our 3.6/im LFs 



are more directly comparable with near-IR LFs at the slightly 
shorter wavelengths of H (1.6/im) and K (2.2//m) as they 
also probe the stellar mass of the galaxies. Figure 12 summa- 



rizes the properties of this and previous near-IR surveys, and 
demonstrates that the IRAC survey covers a previously un- 
explored regime of a large surface area combined with good 
photometric depth. 



The first i - C-band Coma LF was measured by Mobasher 
& Trentham ( 1998 1 for a small 41 arcmin-^ area in the core 



down to = 18.5. This LF was not well constrained 
(a = -lAltoi7 with M^ fixed at -24.35) due to large 
uncertainties in the background source subtraction, but the 
faint galaxies were shown to be dwarf spheroidals by their 
optical-IR colors. A better comparison comes from the larger 
area (660 arcmin^) H-hand study of the core region by de] 
Propris et al. (1998i, where the shape of the LF (down to 



H ~ 16) is similar in many ways to the 3.6/im LF presented 
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TABLE 4 

Other measured slopes of the Coma LF 



Filter 


a 


Magnitude Range" 


Field Size 
(arcmin^) 


Backgroimd Subtraction Method 


Background Field Size 
(arcmin^) 


Reference 
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References. — (1) |Andreon|(T999t; (2)|Cortese et al.||2003); (3)|Beiiersbergen etaLH2002t; (4)|Adami et"arH2007); (5)|Thompson & GregoryHT993); 

(6)'Biviano et al. llggTi; (7) Mobasher et al. (2003); (8) Andreon & Cuillandre 12002); (9) Lobo et af f 1997); (lO)'lglesias-Paramo et al.'f2003^; ( ll)|lglesias-| 
|P|ramo et al. (20021; (12) Trentham ( 19981; (13) Bernstein et al. 1 19951; (14) Lopez-Cruz et al. (19971; (15) Seeker et al. ( 1997); (16) Milne et alT| j2007t ;"(T7) 
[HePropris et al.|^1998^ ; (18)^Andreon & Pello l2000l; (19) ,Mobashe r & Trentham, tl998; ; (20) This work; (21) jBai et al.. (20Q6l. 

Note. — " Apparent magnitude range of a slope measurement. Where only faint data were fitted, or where two-component models were used, the quoted 
magnitude range and slope is for the faint component only (denoted by "F"). Otherwise, slopes were measured over the full LF magnitude range. 
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Fig. 12. — A comparison of previous near-IR coma surveys with the Spitzer 
IRAC survey. All flux densities for the H- and _R'-band surveys were trans- 
formed to equivalent 3.6/im flux densities assuming an elliptical S EP. The 
IRAC su rvey is ~2 times deeper than the deepest //-band survey of Andreonj 
|& Pell6H2000^ , and covers and area 7 times greater than any previous study. 



here. The statis tics (^20 galaxies compared to only ^ 40 
in the [ Mobasher & Trentham 1998 study) allowed better con- 
straints on the parameters, and the bright end was fitted with a 
Schechter function wi th a = -0.78 and = -23.87 (for 
DM = 35.0). Note that |de Propris et al.| ( [T998l ) also found evi- 
dence of an excess of galaxies fainter than H ~ 14.5 mag; us- 
ing two different background subtraction methods (statistical 
field subtraction and B — R color selection), they measured a 



power-law slope of a — —1.7. Another i?-band study by An- 



dreon & Pell6| ( |2000| l covered a smaller area (480 arcmin^) 
15' North -East from the cluster center, and reached on e 
magnitude deeper (H ^ 17) than de Propris et al. ( 1998| l. 



This iJ-band LF was fitted with a single Schechter function 

with a slope of a 1.3 and M;^ = -23.9 (for DM = 

35.0), but again was not well-constrained. A dip was also 
seen at ~ —22, which is also present in both fields in 
the 3.6^m LFs and in the de Propris et al. (1998) i/-band 
data. We also see another shallow dip in the IRAC data in 
both fields at M3 g^m ~ —23.75. These dips have been seen 



in optical surveys of Coma ( Godwin et al.p 983; Biviano et 
[ar][T995j [Mobasher et"aL||2003 '. 'Adami et al. 2007j), and are 
thought to be a result of the merging and disappearance of 



bright elliptical galaxies ( Andreon & Pello [2000 [[Adami et al. 
[20f7] l. 

All of these near-IR LFs have relatively flat faint-end 
slopes, similar to the IRAC 3.6/xm LF when fitted by a 
Schechter function over the range —27 < Ma .g^m < —20.5. 
However, apart from the [de Propris et al. ( [199 8) i/-band 
study, they do not show evidence of the excess population of 
faint (M3 > —20.5) galaxies seen here. This is likely to 
be a direct consequence of the increased sensitivity of IRAC 
(i.e. darker sky), permitting much better statistics compared 
to the ground-based instruments (see § [TJ. The vastly lower 
sky background level at 3.6/im will also increase detection of 
low-surface-brightness (LSB) dwarf galaxies that may have 



been missed in the ground-based surveys. 

Recently, there has also been a longer-wavelength Spitzer 
Multiband Imaging Photometer (MIPS) 24/im study of the 
Coma LF (Bai et al. 2006). However, due to the larger PSF 
of MIPS (6" FWHM), it is much more prone to confusion is- 
sues in crowded fields like Coma and cannot therefore go as 
deep as IRAC. Bai et al. (20061 detected 217 member galax- 
ies at 24/im, but only down to m24^m — 11 (equivalent to 
ni3.6Aim ^15 assuming an elliptical SED). Since the MIPS 
data were not deep enough to reach the dwarf galaxy popu- 
lation turn-up seen in the IRAC data, only single Schechter 
functions were required to fit the data, with slopes of a = 



-0.99^0:52 in the core and a = -1.32^g {4 in their off-center 
region. These results are consistent with the 3.6pm bright-end 
slopes in the two fields, within the combined errors. 

7.3. Comparison with Field LFs 

We can also compare the space density of 3.6/im detected 
galaxies in the Coma LF to that of the fi eld p opulation over the 
same magnitude range. As shown in § 3.3 in addition to the 
obvious excess at bright magnitudes, there is also a significant 
excess of galaxies in Coma compared with 3.6/im field galaxy 
counts of Fazio et al. (2004b) at faint magnitudes (ma g/^m > 
14.5), the position of the Coma upturn (see Figure|7]). 

Recent determinations of field i^-band LFs are deep 
enough to compare with Coma (e.g. IBolzonella et al.||2002| 



1+0.11 



[Pozzetti et a l. 2003 , Feulner et al. 2003 ; Saracco et al. 2006; 
and they are all consistent with a flat faint-end slope for the LF 
(a ^ -1.0 to -1.3) out to redshifts of 2: < 1.3. A 3.6/x m field 
LF has also been published by Babbedge et al. ( 2006) for the 
Spitzer SWIRE survey, which is deeper than our Coma survey 
(m^ Q^jn ^19.5). This also shows a flat faint-end slope, with 
a ^ -0.9 to -1.0 over the redshift range z ~ to 1.5. The 
excess of galaxies detected in Coma with respect to 3.6/im 
field galaxy counts, together with relatively flat slopes of near- 
IR field LFs extending beyond the depth of our IRAC Coma 
sample, support the evidence that the large dwarf population 
detected at 3.6pm is not an artifact of incorrect background 
source subtraction, and that they do indeed belong to Coma. 

7.4. Comparisons with Coma LFs at other wavelengths 

The optical LF studies of Coma are deep enough to reach 
the dwarf galaxy population, and a wide range of faint-end 
slopes have been measured (a ^ -1.0 to -1.8, see Table[4l e.g. 



'Thompson & Gregory' 1993 ; Bernstein et al.'1995VBiviano et 
al. 1995; Lopez -Cruz et al. 1997; Seeker et al. 1997; Beijer" 



sbergen et al. 2002; Andreon & Cuillandre 2002; Lobo et al. 
,1997, .Trentham„1998) . A more recent HST i?-band study 
measures an even steeper slope of -2.3 in a small 5 arcmin^ 
core region ( Milne et al.|2007| l. 



We can qualitati vely compare the 3.6pm LF to the R and B- 
band optical LFs of Mobasher et al.[ ('2003), as the background 
subtraction was also performed using spectroscopic member- 
ship fractions. To translate the optical LFs to their expected 
magnitudes at 3.6pm, we compute mean optical-IR colors of 
the spectroscopically confirmed members to be R-3.6pm ^ 
2.7 and B-3.6pm ^ 4.2 for all galaxies (without regard to 
galaxy classification). However, due to the known luminos- 
ity/galaxy type segregation in clusters, the corrections will not 
be uniform across all magnitude bins. In order to match the 
optical and 3.6pm LFs more accurately, we derive the mean 
R — 3.6/im and B — 3.6/im colors of the confirmed members 
in half-magnitude bins, corresponding to the binning of the 
optical and IR data. 
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Fig. 13. — A qualitative comparison of the combined 3.6^m LF (Coma 1 + Coma 3) with the equivalent R- and B-band optical LFs of |Mobasher et al. ] j2003) . 
These are directly comparable to the IRAC LF as the same method of background source subtraction was used (i.e. spectroscopic membership traction). The 
optical LFs have been shifted to the 3.6/^m frame using the mean optical-lR colors of the confirmed cluster members per half-magnitude bin. The range of stellar 
masses covered by these observations is also shown, assuming a iC-band mass-to-light ratio of 0.91 at z = (Drory et al. 20041. The 3.6^m and optical LFs 
agree well at the bright end up to Ma.e^jm ~ —20.5. At this point, equivalent to > —18.0, there is a steep rise in the IR LF towards fainter magnitudes 
compared to the optical, presumably with the onset of the dwarf galaxy population. Note that while one can compare the shapes of the optical and IRAC LFs, 
a comparison of their relative normalizations may be misleading; the excess of IR-detected galaxies may represent a LSB dwarf population, or a population too 
red to be detected down to the optical limit of i? ~ 20 (or a combination of both). 



Figure 13 shows the combined (Coma 1 + Coma 3) 3.6/im 



LF, together with the equivalent R- and S-band LFs of 
Mobasher et al. ( 2003[ ) shifted by the mean R — 3.6/im and 



B — 3.6fim colors. Also shown is the range of stellar masses 
covered by the IRAC observations, which are estimated as 
M*(Mq) ~ 5.7 X 10^/3. 6^m, where /a.e^m is the 3.6^mflux 
density in /iJy. This is calculated using a distance of lOOMpc 
and a if -band mass-to-light ratio of 0.91 at z = (Drory et al. 
|2004), after transforming the 3.6/im flux densities to Ji'-band 
using the Spitzer EX-PET^ tool, assuming an elliptical SED. 

The bright-ends of the LFs agree well, which is to be ex- 
pected as we are sampling almost the same galaxy popula- 
tions. However, fainter than Ma.g^^m ^^-^-20. 5, the shapes dif- 
fer significantly, with the IRAC LF having a much steeper 
slope than the optical LFs. At the faint end, we show both 

^ http://ssc.spitzercaltech.edu/tools/expet/ 



the raw 3.6//m LF (with no optical selection) and the two 
g — r color-cut LFs, illustrating that the IR-excess is still sig- 
nificant with these constraints. It is, however, important to 
note that while one can compare the shapes of the optical and 
IRAC LFs, a comparison of their relative normalizations may 
be misleading. A possible explanation for the difference be- 
tween the optical and IR LFs is that the apparent excess of 
faint IR-detected galaxies may represent either a LSB dwarf 
population that IRAC is more sensitive to, or a populat ion of 
dwarf galaxies that are too red to be detected in Mobasher et 
al. (2003 ) down to the optical limit of i? ~ 20 (or a combina- 
tion of both). 

However, there are many other optical surveys that probe 
much deeper into the dwarf population, which, by necessity, 
use statistical background source subtraction (e.g. Bernstein 
(i? < 23.5), |Trendiam|1998| (fi < 24) and |Adami 
{R < 24)). The optical color-cuts we impose on 
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Fig. 14. — The R — 3.6/im versus 3.6/im color-magnitude relation for 
the r < 22.5, g — r < 2 color-selected subset of IRAC-detected galaxies 
matched to the SDSS DR5 catalog. A strong red sequence is formed by the 
Coma members aX R— 3.6/jm ~ 2.5 — 3.5, but note that there is an increase 
in redder members towards fainter magnitudes. There is also a general red- 
dening of all galaxies at faint magnitudes. Although a fraction of these will 
be background galaxies, this is consistent with the suggestion that the large 
IRAC-detected faint population in Coma may be redder than the optically- 
selected faint population. 



our 3.6/im faint-end LFs means that we do not include any op- 
tically undetected galaxies. However, the SDSS catalog only 
reaches r ~ 22.5, so in terms of r-band magnitude, this cut 
may be overly-conservative. There could be additional IRAC- 
detected dwarfs that are optically fainter than this, and indeed 
the deeper optical studies already show evidence of such faint 
populations. For example, the i?-band membership fraction of 
[Mobasher etaL] ( |2003| l imphes a flat LF slope out to R 20, 
but the deeper j?-band LFs of 'Bernstein et al. ( 1995 1, Tren- 



[diam, ( ,1998j and |Adami et al.| (2007 ) all continue to rise be- 
yond this with slopes of a ~ -1.4 to -1.7. This may therefore 
be the i?-band regime where the excess 3.6/im dwarf galaxy 
population is detectable, and suggests that this near-IR sur- 
vey may manifest this population. To investigate this further, 
we have cross-correlated the Coma 1 IRAC catalog with the 
deep |Adami et al.| ( |2007| l optical catalog over an area of 1365 
arcmin^ in the central region common to both surveys. We 
find that ^93% of the faint 3.6/im sources (ms e^^m > 14.5) 
have optical matches down to i? = 24. 

These results imply that the faint IR population has much 
redder optical/IR colors compared to the optically selected 
dwarfs in the cluster Since we have used membership frac- 
tions to construct the LF, we have not yet identified which par- 
ticular galaxies in this population are Coma members. How- 
ever, in Figure 14 we show the r — 3.6/im versus 3.6/im color- 
magnitude relation for the r < 22.5, g — r < 2 optically- 
selected subset of IRAC-detected galaxies matched to the 
SDSS catalog. In a similar fashion to the optical color dia- 
the Coma members form a strong red se- 
.6/im ~ 2.5 — 3.5, but we also note an increase 



gram (Figure [TT i 
quence at r - 

in redder members towards fainter magnitudes. There is also 
a general reddening of all galaxies at faint magnitudes; al- 
though a fraction of these will be background galaxies, this 
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Fig. 15. — Measured values of the Coma LF faint-end slope a from the lit- 
erature as a function of wavelength (see Table|4j. The symbols denote differ- 
ent methods of background source subtraction: control field {black circles), 
B-R color (red squares) and redshifts/membership fractions (blue triangles). 
The 3.6/im 'raw' faint-end slope (with no optical color-cut) is highlighted 
with a large blue open triangle, with an error bar representing the full range 
of a2 obtained with all g — r color cuts in both Coma fields. No signifi- 
cant systematic variation of slope with wavelength is evident in these com- 
bined r esults. The 3. 6/im slope is steeper than all but one measurem ent in the 
UV ( Andreonl 1999k and the HST i?-band core measurement of ( Milne et al.l 
|2007). 

is consistent with the suggestion that the IRAC-detected faint 
population in Coma is redder than the optically-selected dwarf 
population. This will be explored in more detail in a future 
paper. 

7.5. Comparison of faint-end slope versus wavelength 



Figure 15 shows the multi- wavelength Coma faint-end LF 
slopes from Table |4] plotted versus filter wavelength. Al- 
though the areas and photometric depths covered by these 
studies vary, they can be used to investigate any broad trends 
in changes of the faint-end slope with wavelength. Where 
slopes were measured in different regions of Coma in a par- 
ticular study (i.e. core/off-center/total), only the slope of 
the total area LF is plotted. Some evidence has previously 
been found for change of the faint slope between different 
bands (e.g. [Beijersbergen et aL]|2002[ |Adami et al.||2007 



but no significant trend is evident in this figure. The dif- 
ferent symbols/colors denote the different methods of back- 
ground source subtraction used. At optical wavelengths, the 
slopes based on samples measured with a spectroscopic back- 
ground subtraction appear systematically shallower (in agree- 
ment with e.g. Valotto et al.|20 01"Adami et al.'2000'). How- 
ever, the UV measurement of Andreon (1999) and Cortese] 
et al.| (|2003 ), both using a similar method, are much steeper 
(a = -2.16 and -1.65 respectively). The 24/im slope mea- 
sured by |Bai et alT ( 2006 1 using redshift membership frac- 
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tions (a — —1.49) is consistent within the errors with the op- 
tical measurements using statistical background source sub- 
traction. However, this is relatively shallower and therefore 
not directly comparable. A systematic study using reliable 
redshift background source subtraction across all wavebands 
is required to properly determine if there is any color evolu- 
tion of the faint-end slope. Note that the IRAC 3.6/xm slope is 
significantly steeper than all but one measurement i n the UV 
dAndreon 1999) and the deep HST i?-band study of [Milne et 
m H2007l l. 

7.6. Environmental Dependence 

One of the goals of this study is to search for changes in 
the shape of the faint end of the LF in different regions of 
Coma. The shape of the LF is not thought to be universal, 
but strongly influenced b y local environment. For instance, 
Lopez-Cruz et al. ( 1997 1 measure flat slopes {a ^ —1.0) in 



a sample of rich clusters, and find much steeper slopes in 
poorer clusters (—2.0 < a < —1.4). In the central regions 
of clusters, processe s such as tid al stripping and ram pres- 
sure stripping (Gunn & Gott 1972| l conspire to destroy smaller 
dwarf galaxies by removing their stars and gas. Another de- 
structive process that depend s on cluster density is galaxy ha- 
rassment ( jMoore et al.|1998] l, where frequent and rapid close 
encounters of small spiral galaxies convert them into dwarf 
spheroidal systems. In addition, Adami et al. ((2000) sug- 



gest that the high temperatures present in Coma during its ini- 
tial formation could have suppressed the formation of dwarf 
galaxies in the core region. These effects point towards a sce- 
nario where the number of dwarf galaxies should decrease 
towards the center of the cluster, leading to a flattening of 
the LF faint-end slope. However there are other competing 
factors that come into play. The denser the environment, the 
more likely that dark halos of galaxies could collect gas from 
the inter-cluster medium to be converted into stars, leading to 
larger numbers of dwarf galaxies in denser evolved cores than 
the le ss-dense, unevolved off-center or field (Mobashe r et al.| 
|2003| and references therein). This could counter-balance the 
effects of destructive forces of tidal interactions, ram pressure 
stripping and harassment. 

We do find significant differences between slopes at the 
core and off-center regions of the cluster If this directly traces 
the mass, then a difference in the LFs means galaxies have 
different formation/destruction and merger histories depend- 
ing on their local environment. At the bright end (Schechter 
form), the slope increases from ai = -1.18±0.06 in the core 
to -1.30±0.08 in the off-center field, while at the faint end 
(power-law form) there is an increase from a2 = -2.18±0.13 
to -2.60±0.33 (see Table [3] and Figure [T0|." Additionally, 
Mg is fainter in Coma 3, which reflects the fact that this 
off-center region does not possess as many of the very bright 
galaxies seen in the central region. Some previous studies 
covering sufficiently large areas of Coma have also found 
similar changes in LF slope with environment, a s shown in 
Table [4] For example, Beijersbergen et al. ( 2002| l finds a U- 
band slope of -1.32 in the core region, but a steeper slope of 
-1.54 when all areas of the cluster within the survey are in- 
cluded. A smaller increase is also found in the r-band, with 
no such difference found in the i3-band. A similar effect is 
found in the y-band by |Lobo et al. ( 1997| l. More recently, 
|Adami et al.| ( |2007] l have found significant differences in the 
faint-end slopes of northern and southern areas of the center 
of the cluster, with a steep i?-band slope in the north (-1.74) 
and a flatter slope in the south (-1.28). Mobasher et al. (2003 i 



also find a small increase in the i?-band s lope between the 
core (-1.17) and the off-center (-1.29), and |Bai et all ( |2006l l 
found a similar effect at mid-lR wavelengths, with a slope of 
-0.99 in the core, increasing to -1 .32 in their off-center region. 
Collectively, these results provide strong evidence that dwarf 
galaxies are being destroyed or merged into larger galaxies 
in the denser central regions of the cluster. This scenario is 
further supported by evidence of disruption in the central re- 
gions of Coma in the form of tidal features, resulting from 
galaxy-galaxy and galaxy-cluster interactions (Trentham & 
pobasher. 1 998, .Gregg & West| 1998| ). 



7.7. Redshift evolution of the IR LF 

It is interesting at this stage to compare the Coma IR LF 
with examples at higher redshifts to study possible evolution- 
ary trends. In de Propris et al. ( 1999) , the i^-band LFs of 38 
clusters with redshifts in the range 0.1 < z < 1 were mea- 
sured. Although they were unable to constrain the faint-end 
slopes of the clusters with these data, they found that the evo- 
lution of was consistent with a model where the giant 
early-type galaxies in the clusters formed all their stars in a 
single burst at high redshift (z > 2), and evolved passively 
thereafter. This result was independent of the richness of the 
clusters. 

More recently, there have been two studi es of the IRAC 
3.6/im redshift evolution of clusters. Muzzin et al.| ( 2005] l 
measured the 3.6/xm LFs of 123 galaxy clusters in the Spitzer 
First Look Survey (FLS) in the redshift range 0.15 < z < 
1.22. They found that Mg g evolved in a manner consistent 



with a passive evolution model (in agreement with de Propris 
et al.|[T999) . Moreover, they also found some evidence that 
the faint-end slope became shallower at z > 0.5. Together, 
these results imply similar redshifts (z > 2) for the forma- 
tion of massive galaxies and clusters, with many of the dwarf 
galaxies accreted from the field at z < 0.5. This result is 
also supported by Toft et al. ( 2004| l, whose Kg-hand LF of 
a cluster at z = 1.237 has a faint-end slope of a ~ —0.64, 
which is significantly shallower than clusters in the local Uni- 
verse. Another study by |Andre on (2006) of 32 clusters in the 
Spitzer SWIRE survey in the redshift range 0.2 < z < 1.25 
also found that their values of Mg g^„j were consistent with 
a passive evolution model, but were unable to determine any 
change of slope with redshift with these data. 

The large population of dwarf galaxies detected at 3.6/im 
in Coma supports the scenario of accretion of dwarf galaxies 
from the field at low redshift. However, our values of Mg g^„ 
(-25.17±0.10 for Coma 1 and -24.69±0.15 for Coma 3) are 
slightly brighter than the value predicted by the passive evo- 
lution model of |Bruzual A. & Chariot] ( fTW3l l of -24.33 (A. 
Muzzin, private communication). 

8. SUMMARY & CONCLUSIONS 

We have conducted a wide-field Spitzer IRAC 3.6/im sur- 
vey of two fields in the Coma cluster, in which we have dra- 
matically improved constraints on the faint end of the near-IR 
Coma LF compared to previous IR surveys. Thanks to the ef- 
ficiency of the IRAC instrument and the significantly reduced 
sky background in space, we have been able to probe as faint 
as mg.gpm ^ 17 (Mg.g^m ~ —18) ovcr a large fraction of 
the cluster (~ 1.3 deg^), giving a very representative view 
of the IR dwarf population of Coma. The first of our two 
fields covers an area of 0.733 deg'^ in the core of the cluster 
(Coma 1), and the second is an off-center 0.555 deg^ region, 
-■57' from the core (Coma 3). We detect a total of 29,208 



Spitzer IRAC LF of the Coma Cluster 



17 



sources, 17,872 in Coma 1 and 11,336 in Coma 3. Our main 
results can be summarized as follows: 

1. We construct 3.6/im LFs for the Coma 1 and Coma 3 
fields. To subtract the contribution to the 3.6/im counts 
due to background galaxies, we use spectroscopic red- 
shifts for IRAC-detected galaxies to construct spectro- 
scopic membership fractions for each field. At the 
bright end (M3 g^m < —20.5), the LFs can be fitted 
with a classic Schechter function with slopes of ai = 
-1.18±0.06 (Coma 1) and -1.30±0.08 (Coma 3), and 
M* g^,„ values of -25.17±0.10 and - 24.69±0.15 for 
Coma 1 and 3 respectively. 



2. Faintward of Ma.gpm ~ —20.5, we detect a steep 
faint-end slope in both fields. To guard against any 
optical selection bias in the spectroscopic complete- 
ness functions, we use SDSS photometric data to con- 
struct optically-selected 3.6/im faint-end slopes for 
both fields. Taking only the IRAC sources with optical 
counterparts at r < 22.5 mag and cut at 5 — r < 1.3, 
we measure steep faint-end slopes of 0:2 = -2.18±0.13 
in Coma 1 and -2.60±0.33 in Coma 3, which infer that 
there are at least ~ 1600 member galaxies over the two 
fields. The differences in slopes between the two fields 
provides evidence for environmental effects. 

3. When qualitatively compared to B- and i?-band opti- 
cal LFs constructed from the same areas of Coma and 
using the same spectroscopic completeness method, we 
find that the optical LFs have fairly flat distributions at 
the faint end {a ~ —1.2). A possible explanation for 
the difference between the optical and IR LFs is that 
the excess of faint IR-detected galaxies may represent 
either a LSB dwarf population that IRAC is more sensi- 
tive to, or a population of dwarf galaxies that are too red 
to be detected in the optical survey down to the limit of 
i? 20 (or a combination of both). 

This is the first Spitzer IRAC 3.6/im study of a low-redshift 
cluster, and it has demonstrated that IRACs combination of 
sensitivity, spatial resolution and mapping capabilities makes 
it the current instrument of choice to probe the stellar mass 
distribution of clusters. The steep faint-end rise seen at 3.6/im 
has huge implications towards constraining the mass function 



of galaxies, and provide clues for an optically faint popula- 
tion of low mass galaxies in clusters. Approximately 60% of 
the IRAC sources have optical counterparts in the SDSS data 
down to r 22.5, and > 90% of a subset in the Coma 1 field 
have counterparts down to i? = 24. There is therefore great 
potential that the numbers of IRAC detected galaxies implied 
by our SDSS optical color-cuts as belonging to Coma will in- 
crease with further investigation. 

The chief source of uncertainty at faint magnitudes in stud- 
ies such as this is spectroscopic completeness, which for 
Coma is < 10% at r = 21. To overcome this, we are con- 
ducting wide-field spectroscopic programs to dramatically in- 
crease the coverage of faint sources, which will give us a 
much better constraint on the contribution to the background 
from field galaxies and background clusters/groups. However, 
to cover thousands of objects over square degree fields of view 
(as is needed in Coma) is still observationally very difficult. 
The current state-of-the-art in wide-field multi-object spec- 
troscopy reaches r ^ 22 in several hour exposures (e.g. the 
Hectospec work of Papovich et al."2006), and so at present, 
statistical field subtraction must still be used fainter than this. 

We are also conducting a multi-wavelength study of the 
faint galaxies detected by IRAC, to measure their photometric 
redshifts (and hence their Coma membership), and also to as- 
certain properties such as stellar ages, star formation histories 
and morphological types. Other studies using new GALEX 
UV and XMM-Newton X-ray data of these Coma fields are 
also in progress. Through such multi-wavelength character- 
izations of the galaxy populations of Coma, we will gain a 
better understanding of galaxy evolution in cluster environ- 
ments. 
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